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Abstract 

We calculate, numerically, the low-lying spectrum of closed confining flux tubes that carry 
flux in different representations of SU(iV). We do so for SU(6) at (3 = 171, where the calculated 
low-energy physics is very close to the continuum limit and, in many respects, also close to 
N = oo. We focus on the adjoint, 84, 120 , k = 2A,2S and k = 3A, 3M, 3S representations 
and provide evidence that the corresponding flux tubes, albeit mostly unstable, do in fact 
exist. We observe that the ground state of a flux tube with momentum along its axis appears 
to be well defined in all cases and is well described by the Nambu-Goto free string spectrum, 
all the way down to very small lengths, just as it is for flux tubes carrying fundamental flux. 
Excited states, however, typically show very much larger deviations from Nambu-Goto than 
the corresponding excitations of fundamental flux tubes and, indeed, cannot be extracted in 
many cases. We discuss whether what we are seeing here are separate stringy and massive 
modes or simply large corrections to energy levels that will become string-like at larger lengths. 
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1 Introduction 

In the confining phase of SU(iV) gauge theories in 3 + 1 or 2 + 1 dimensions, the flux between 
sources in the fundamental representation is carried by a flux tube that at large separations, 
I, will look like a thin string. The spectrum of such a string-like flux tube, whether closed 
(around a spatial torus) or open (ending at two sources), should be calculable from an effective 
string action [fl [2] once I is large enough that the energy gap to the ground state has become 
small compared to the gauge theory's dynamical scale, ~ O(Aj-fg) in D = 3 + 1 and ~ 0(g 2 ) 
in D — 2 + 1. Indeed, it may be that the spectrum is simple even at smaller /, where the 
energy gaps are large, once N is so large that flux tubes effectively do not mix or decay. 
In recent years a great deal of progress has been made in determining the universal terms 
of this effective string action thus determining the spectrum at large /. (See [3] for a recent 
review.) Simultaneously, numerical lattice calculations have determined the spectrum at small 
to medium values of /, where the dynamics turns out to be remarkably close to that of a free 
string theory (Nambu-Goto in flat space-time). 

In this paper we extend our recent lattice calculations of the spectrum of closed flux 
tubes in 2+1 dimensions [6] to the case where the flux is in representations other than the 
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fundamental. This will include cases where the flux tube is stable for all I (e.g. the ground 
states of A/"-ality k = 2 or k = 3) and cases where it is not. Whether the latter do have a well- 
defined identity is an interesting question which we shall also address, albeit only empirically 
in this paper. (We are not aware of a quantitative theoretical analysis of the binding and decay 
of such 'composite' flux tubes, although the general framework for decays has been developed 
in [1], and it would be interesting to understand if the flux tubes considered in this paper 
satisfy the conditions for those calculations to be accurate. See also [5] for related work.) As 
in our earlier work [6] nearly all our calculations are in SU(6), where the theory is close to 
its N = oo limit for many low-energy quantities, but far enough away from that limit for the 
k = 2 and k = 3 ground states to be well below their decay thresholds. Our calculations are 
at a fixed value of the lattice spacing a that is small enough for most lattice corrections to be 
negligible (within our statistical accuracy). 

In the next Section we provide a (very) brief sketch of relevant analytic and numerical 
results. We then describe the technical aspects of the lattice calculation. In Section @] we 
present our results. We begin with flux tubes carrying flux in the k = 2 symmetric and anti- 
symmetric representations (that arise from /®/, where / is the fundamental representation), 
then move on to the three minimal k = 3 representations (arising from /<g) /(g) /), the adjoint 
flux tube (from / <g> /) and those carrying flux in the 84 and 120 representations (arising from 
/ <S> f <S> /). The Appendix describes the properties of these representations. Such flux tubes, 
when they exist, can be thought of as bound states of (ant i) fundamental flux tubes and their 
spectra should contain the imprint of the massive modes associated with that binding. The 
latter should be additional to the usual massless stringy modes, which are the only ones to 
appear in the spectrum of fundamental flux tubes in D = 2 + 1 [6 J . 

The lattice calculations are very similar to our earlier work with fundamental flux and we 
refer to that work [6] for most of the technical details. We also note our earlier calculation 
of the spectrum of k = 2 flux tubes [7] performed at smaller N and for coarser a, and to 
earlier calculations of fc-string tensions |8|. We refer to these for a more detailed discussion of 
fc-strings. 

2 Background and Overview 

We are interested in the spectrum of flux tubes that are closed around a spatial torus of length 
I. We make the sizes of the transverse spatial torus, l±, and the (Euclidean) temporal torus, l t , 
large enough that the resulting finite size corrections are negligible. As I decreases, the theory 
suffers a finite volume transition at I — l c — 1/T C where T c is the deconfining temperature, 
and for I < l c the theory does not support winding flux tubes. This transition is strongly first 
order for SU(6), the case of interest in this paper. Since T c ~ y/af in terms of the fundamental 
string tension, this means we can study closed flux tubes of length I > l/-y/<7/. 

Since the spectrum of the Nambu-Goto model turns out to be an excellent starting point 
for much of the observed fundamental flux tube spectrum, we begin by briefly summarising 
it. We then say something about relevant analytic results for long flux tubes - an area in 
which striking progress has been made in the last few years - as well as the numerical results 
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for flux tubes in the fundamental representation, which the present work extends to higher 
representations. We then say something about those higher representations. 

2.1 analytic expectations 

Recall that we consider flux tubes that are closed around a spatial torus of length I, with 
the transverse and Euclidean time tori chosen so large as to be effectively infinite. Such flux 
tubes may carry non-zero longitudinal momentum. (We do not consider non-zero momentum 
transverse to the string since that does not teach us anything new.) In the iV — > oo limit where 
decays and mixings are suppressed, the world sheet swept out by the propagating flux tube has 
no handles or branchings and so has the simple topology of a cylinder. The simplest effective 
string action is proportional to the invariant area of the sheet in flat space-time (Nambu- 
Goto). The Nambu-Goto spectrum arises from left and right moving massless 'phonons' on 
the background string of tension a. Let n^Ft){k) be the number of left(right) moving phonons 
of momentum \p\ = 2ftk/l and define their total energy to be 2ttNl(r)/1, i.e. 

N L = Y, n L(k)k, N R = Y,n R {k)k. (I) 

k k 

so that the state has total longitudinal momentum p = 2nq/l with 

N L -N R = q. (2) 
The energy levels in D = 2 + 1 turn out to be given by P, \TU\ 

WH{(^ + 8™(^4) + (^) 2 }' (3) 

where n = Nl + Nr and one can readily calculate the degeneracy of a given energy level. We 
note that the parity of a state is given by 

p ^ -^number of phonons (4) 

We display in Tabled] the states which we will later discuss in more detail. (Here a±k creates 
a phonon of momentum ±2irk/l.) We refer to [6] for a more detailed discussion, and reasons 
why we ignore quantum numbers other than parity and momentum along the flux tube. 

Note that the spectrum in eqn([3]) is derived using naive light-cone quantisation |9J. Its 
actual relationship with Nambu-Goto in D < 26 is a subtle question, which is considered 
critically in We will nonetheless use it for comparative purposes and refer to it as 'the 
Nambu-Goto spectrum'. 

For large enough I we can expand eqnQ in powers of 1/al 2 . The first correction to the 
linear al piece coincides with the well-known 0(1/1) universal Liischer correction [TJ [2]. It is 
now known that the 0(1/Z 3 ) correction is also universal [TO} ITT] as is the 0(l/l 5 ) correction 
|12] . (This is for D = 2 + 1; there are interesting differences in D = 3 + 1 [3]). The universality 
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Table 1: The states of the lowest Nambu-Goto energy levels with p = 2nq/l for q = 0, 1,2, 
and g = excited states with N L + iV R < 6. 

class is determined by the massless modes living on the string. If, as is plausible here, the 
only such modes are those arising from the bosonic massless transverse oscillations, then these 
universal terms coincide with the corresponding terms in the expansion of the Nambu-Goto 
action and energy levels [T2], [13]. Thus once / is large enough for the expansion of eqnQ in 
powers of l/l 2 a to converge (which occurs at small / only for the absolute ground state) we 
can expect the free string Nambu-Goto theory to provide an increasingly accurate description 
of that part of the closed flux tube spectrum. 

Note that such effective string calculations become valid for an excited flux tube once / 
becomes large enough that the energy gap AE n ~ E n (l) — ol becomes small compared to the 
dynamical energy scale of the theory ~ -y/c. And this is so independent of N. However the 
expansion of eqn([3]) only requires AE n < ol which is a weaker condition. So for the effective 
string approach to be valid all the way down to the Nambu-Goto radius of convergence we 
presumably need to invoke nearness to the N = oo limit as well. 

While the above analytic progress has so far concerned flux tubes at large enough I, we 
remark that there have been promising recent attempts at understanding the spectrum at 
smaller / from considerations of the scattering matrix of phonons on the world sheet [H] . (See 
also [15J.) 

These analytic results assume that the flux is carried by a single flux tube. While this 
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is indeed the case for fundamental flux tubes, in appropriate limits, it is not clear what 
happens for higher representations R. While we can still expect an effective (Goldstone) 
action approach to be valid as long as AE n ~ E n (l) — ctrI < y/crf, extending the range of 
validity by an appeal to large N is dubious. Indeed, in the N = oo limit we expect the flux to 
be carried by an appropriate number of non-interacting fundamental flux tubes. (As we shall 
see below when we consider explicit operators for such flux.) Thus we are not able to rely on 
an ideal N — > oo limit in the same way as we can for fundamental flux tubes. 

2.2 fundamental flux tubes 

In [6] we performed calculations in SU(6) of the closed flux tube spectrum on the same lattices, 
and at the same coupling as in this paper. We briefly list some of the conclusions of that work 
that are relevant to this paper. 

1) The absolute ground state is very accurately described by the free string prediction in 
eqnQ, with a correction only becoming visible for l\/(Jj < 2. 

2) This correction is consistent, within the errors, with being either oc l/l 5 or oc l/l 7 , where 
the latter is the prediction of the analysis of universal terms. 

3) The lightest states with p ^ also show no visible correction to Nambu-Goto down to 
ly/(jf ~ 1.5. These states contain phonon excitations and so we see that the flux tube behaves 
like an excited thin string even when its length is about the same as its width (which is naively 

4) In general whenever an excited state corresponds to phonons that are all right or left mov- 
ing, corrections to Nambu-Goto are almost invisible. 

5) While other low-lying excited states typically show larger corrections, these typically be- 
come insignificant at values of I that are much smaller than required for the expansion of 
eqnQ in powers of l/l 2 a to become convergent. That is to say, our results show that the 
Nambu-Goto prediction is still good when all the terms in the l/l 2 a expansion are important: 
i.e. the series of correction terms must itself resum to a modest total correction even at small 
I. 

6) There is no evidence at all of any non-stringy massive modes that are additional to the 
stringy ones that are well described by the free string theory spectrum. 

One of our main motivations for the present study is to contrast the above with what one 
finds for flux tubes that are bound states of fundamental flux tubes, and where the binding, 
measurable through the value of the string tension, provides unambiguous massive dynamics 
that should somehow make itself seen in the flux tube spectrum. 

2.3 flux tubes in higher representations 

Consider two well separated sources in representations R and R. The flux between them will 
be carried by one or more confining flux tubes and, if we ignore the possibility of screening, 
will be in the representation R. 

The representations of SU(6) that we consider in this paper are the fundamental /, the 
adjoint A which appears in / <g> f, the representations 84 and 120 which appear in / <g) / <g) /, 
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and the various irreducible representations generated by / <8> / and / <g) / (g) /. These last 
two belong to the k = 2 and = 3 sectors respectively. That is to say under a global 
gauge transformation that is an element of the centre, g(x) = e m ^ N I, the sources transform 
as = e mk l N . Under this categorisation the /, 84 and 120 belong to k = 1 and the adjoint 
A to k = 0. In the k = 2 sector we consider the antisymmetric 2A and symmetric 25* 
representations. In the k = 3 sector we consider the antisymmetric 3A, the mixed 3M, and 
the symmetric 3S representations. All these representations are discussed in more detail in 
the Appendix. 

Since gluons transform trivially under the centre, screening cannot change the value of k. 
Hence the absolute ground state in each fc-sector will correspond to an absolutely stable flux 
tube. These are often referred to as fc-strings, although this term is often used more loosely 
to label all states in a given fc-sector. Note that there will be an absolutely stable ground 
state for each parity, P, and longitudinal momentum, p, within each /c-sector. (Note also that 
at a given / the lightest state with such non-trivial quantum numbers may include a glueball 
that carries some of the quantum numbers. Such states decouple from our calculations in the 
iV — > oo limit and, as we shall see, appear to play no role even for N = 6.) 

Earlier work [161 El U\ has shown that the k = 2, 3 ground states are almost exactly 2 A 
and 3A respectively, except when the flux tube is very short, I ~ l c . This is related to the 
observation that, despite the fact that gluon screening can take one from e.g. 2S to 2A, the 
actual overlap is found to be extremely small [8J. This interesting feature of the dynamics is 
something we shall examine in more detail in this paper. 

We note that some overlaps are lower order 'ml/N and hence would be naturally suppressed 
for SU(6). This includes the overlap of the adjoint flux tube to the vacuum (or glueballs) and 
the 84 and 120 flux tubes to a single fundamental, /, flux tube. On the other hand the overlap 
of the adjoint onto a pair of flux tubes, one / and the other /, should not be suppressed. 
Similarly for 84 and 120 to three flux tubes, 2 fs and one /. We will be careful to discuss 
these possibilities when we present our results below. 

3 Lattice methods 
3.1 lattice setup 

Our space-time is a periodic cubic L x x L y x L t lattice with lattice spacing a. The degrees of 
freedom are SU(iV) matrices, U^(x,y,t) or more compactly Ui, assigned to the links I of the 
lattice. The action is the standard Wilson plaquette action, so the partition function is 




(5) 



where U p is the ordered product of matrices around the boundary of the elementary square 
(plaquette) labelled by p. Taking the continuum limit, one finds that 



where g 2 is the coupling and ag 2 is the dimensionless coupling on the length scale a. The 
continuum limit is approached by tuning /3 = IN jag 2 — > oo. 



3.2 calculating energies 

Here we give a brief sketch and refer the reader to Section 3 of [6] for a detailed exposition. 
We calculate energies from the time behaviour of correlators of suitable operators {4>i}, 

c t3 {t) = (4mm = (<py Hant <p 3 ) = E c ^~ aEknt - (7) 

k 

Since we wish to project onto loops of flux closed around the x-torus, we use operators that 
wind around the x-torus. The simplest such operator is the Polyakov loop 

l p {n y ,rit) = [[ U x (n x ,n y ,nt) ; (p(n y , n t ) = Tr R {l p (n y , n t )} (8) 

n x =l 

where I = aL x (we shall measure I in physical units and L in lattice units, unless indicated 
otherwise) and we have taken the product of the link matrices in the x-direction, around the 
x-torus and the trace is taken in the desired representation R. We also use many other winding 
paths, as listed in Table 2 of [6], and also with smeared and blocked SU(iV) link matrices [6]. 
Using all these paths we can project onto different longitudinal momenta and parities. The 
transverse momentum dependence is determined by Lorentz invariance and so we only consider 
p± = operators, obtained by summing over spatial sites, e.g. l p (pi_ = 0,n t ) oc Ylm y lp( n yi n t) 
in eqn([H]). Unless otherwise stated all winding operators in this paper will be with p± = 0. 

We now perform a variational calculation of the spectrum, maximising (e~ Ht ) over this 
basis (usually projected onto the desired quantum numbers) . We usually do so for t = a and 
this provides us with an ordered set of approximate energy eigenoperators {ipi}- We then form 
the correlators of these, (^(^^(O)), and extract the energies from plateaux in the effective 
energies, defined by 

<WM*(Q)> =exp{ - agi , //(t)} . (9) 

These plateaux typically begin at values of t that are larger than t = a. Given the propagation 
of statistical errors, we can only identify such a plateau if it corresponds to the operator having 
a large overlap onto the desired state. Note that this largely excludes the possibility that our 
energy estimate is contaminated by a small admixture of a lower lying state. (The effective 
energy only provides an upper bound on the desired energy if extracted where we perform the 
variational calculation, i.e. t = a in our case.) It is only where we have significant evidence 
for a plateau that we quote an energy. 

This above procedure is appropriate for stable states. However many of our states will 
be unstable. (We will usually indicate that in our figures.) If these states are analogous to 
narrow resonances then they are just as relevant to us as they would be if stable. If the 
decay width is very small (as it often might be because iV is quite large) then by continuity 
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Table 2: String tensions for various representations (see text) in SU(6) at = 171. Also ratios 
to the fundamental and predictions of Casimir scaling. 



we expect that within our finite errors the correlators will behave just as they do for stable 
states. Conversely, if our correlator looks just like that of a stable state, with an apparently 
well-defined energy plateau, we can assume that the state is very narrow, and extract an 
energy. This will certainly not always be the case. Sometimes we have accurate correlators 
out to large n t where there is no sign of a plateau, presumably because the state has a large 
decay width. We shall perform a heuristic analysis of some of these cases when we come to 
them. The interesting conclusion will be that this leads to an energy much higher than one 
would naively guess by looking at the effective energies. 

We remark that the exact eigenstates of H consist of asymptotic states composed of any 
stable flux tubes and scattering states of these. (And in addition, at finite N, of stable glue- 
balls.) In particular this includes scattering states of fundamental and antifundamental flux 
tubes with various relative momenta. However our basis of operators will usually (although 
not intentionally) have a small overlap on these, and so we usually will not see them in our cal- 
culation. We will comment further on this when we consider examples of what are presumably 
unstable states. 

4 Spectrum results 

In this section we present our results. Before entering into details we list in Table [2] the 
string tension gr that we obtain by fitting the absolute ground state energy, E Q (l), for each 
representation with the Nambu-Goto expression in eqn ffTOl plus a 0(1// 7 ) correction, i.e. 

i 

^w^'l 1 -!^}^ < 10 > 

We use this correction because it is the leading correction to the universal terms j3] , but since 
any correction will only affect Eq(1) at small I our particular choice does not affect the value 
of the extracted string tension. We compare the ratio gr/gj to the ratio of the quadratic 
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Table 3: The energies, E(l), of the lightest flux tubes of length / (all I in lattice units) and 
p = 0, and with the flux belonging to the indicated representations. The fundamental (f) 
values are from [Sj. For SU(6) at j3 = 171.0. 



Casimirs, Cn/Cf. There are old arguments for such 'Casimir scaling' (see [T6J for a discussion 
and references) as well as newer ones, e.g. |18j . We see from Table [2] that it works remarkably 
well. This corroborates earlier studies [T7"t |8] for some of these representations. (As well as 
older studies in SU(3) of open flux tubes, e.g. [19].) The values of E (l) that go into these 
fits are listed in Tables [3] and 0] where we also show the lattice sizes used. For completeness 
we include the values for the fundamental representation obtained in our earlier work [6]. 

4.1 finite volume corrections 

Calculations on I x lj_ x l t lattices will suffer finite volume corrections if l± and l t are not large 
enough. This problem becomes more severe as I decreases. Some checks have been performed 
in [211 [7] for = 2 flux tubes, and in [6] for excited states of k — 1 flux tubes. Since our 
calculations are now more accurate, it is worth revisiting this question. 

We focus on our shortest flux tube, where we employ a 16 x 100 x 200 lattice. We are 
confident that l t /a — 200 is long enough since e~ Elt = e - 200aE ( l=ie ) [ s negligible for all the 
l/a = 16 flux tube energies listed in Tables [3] and |U We therefore test whether l±/a = 100 
is large enough and we do this by performing calculations on 16 x l±/a x 200 lattices with 
l±/a = 20,40,60,80. To speed up these very slow calculations we use a much reduced basis 
of operators - just the simplest Polyakov loops at various blocking levels. This still allows us 
to obtain accurate values for the ground states but not for any of the excited states. (Which 
is why we introduced our extended operator basis in the first place.) So for the excited states 
we continue to rely on the study in [6] and the rescaling of those results to our lattice spacing. 

In Table [5] we show our results for the ground states in various representations. We see 
that the fundamental flux tube suffers no finite volume corrections for lj_/a > 40 within the 
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0.954(38) 



Table 4: The energies, E(l), of the lightest flux tubes of length I (all I in lattice units) and 
p — 0, and with the flux belonging to the indicated representations. For SU(6) at (3 = 171.0. 



statistical uncertainty of about ±1%. For the higher representations there are still visible 
corrections for lj_/a = 40 but lj_/a > 60 appears to be safe at the ±2 or 3 percent level of our 
statistical errors. It thus appears that l±/a = 100 is in fact a very safe and conservative choice. 
This provides further evidence that the energies calculated in this paper are not afflicted by 
significant finite size corrections. 

In Table M we again show some results for the ground states (and also for some excited 
states) in various representations, but this time for a much longer flux tube, l/a = 40. This 
confirms that a transverse size l±/a = 32 is already large enough, and the sizes we have 
actually used are very conservative. 

Finite size corrections also affect the screening of one representation to another, as shown 
in Tables 2,3 of [7]. This is relevant because it is only when the screening is very weak that 
we can categorise the states as being (almost entirely) in k = 2A and k = 2S rather than just 
k = 2 (and similarly for our other representations). We therefore perform a similar analysis 
here. We define the normalised overlap 

n , M (gUg = 0) ^ b{t = 0)) 

OoAS (0) = J ■ 1 11 

<«W* = 0)W* = W 2 (*lsj,(t = 0)*2s fi (t = 0)) 1/2 

where <&r(£) is the simple Polyakov loop at blocking level b and representation R, in the time- 
slice t — and, as usual, summed over spatial sites so as to have zero transverse momentum. 
(Obviously we will average over all equal times.) The range of values of b is restricted by the 
fact that a 'blocked link' [221 E] joins lattice sites that are separated by 2 b_1 lattice sites. So 
for l/a = 16 it only makes sense to consider 1 < b < 5. Essentially, loops at blocking level b 
are smeared over distances significantly greater than this separation 2 6 " 1 . Thus the highest 
blocking level shown typically involves operators that overlap over the boundary of the torus 
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aE R (l = 


16; l±) 






R 


l± = 20 


l± = 40 


l± = 60 


l± = 80 


l± = 100 


f (k=l) 


0.0742(10) 


0.0781(8) 


0.0781(11) 


0.0766(13) 


0.0777(3) 


k=2A 


0.1167(18) 


0.1385(18) 


0.1430(28) 


0.1460(20) 


0.1460(14) 


k=2S 


0.2335(18) 


0.2243(19) 


0.2260(31) 


0.2280(24) 


0.2256(29) 


k=3A 


0.1292(32) 


0.1624(26) 


0.1706(35) 


0.1748(26) 


0.1742(11) 


k=3M 


0.2521(48) 


0.2573(40) 


0.2675(64) 


0.2699(60) 


0.2705(61) 


k=3S 


0.4148(41) 


0.364(12) 


0.390(9) 


0.421(5) 


0.409(4) 


Adj (k=0) 


0.1553(44) 


0.1652(30) 


0.1692(47) 


0.1796(44) 


0.1658(66) 



Table 5: The energy, Er(1), of the lightest flux tube of length I = 16 (all I in lattice units) on 
a 16 x /j_ x 200 lattices, and with the flux belonging to the representation R. For SU(6) at 
P = 171.0. 



and these can be affected by strong finite volume corrections. 

Bearing the above in mind, we show our results for the overlap 02A,s{b) in Table [71 We 
remark that the calculations with l/a ^ 16 are mostly with lower statistics, designed to be 
sufficient for our purposes here. We also calculate Polyakov loops in the (usually) longer y 
direction, and this gives us some values of C^s^fe) at small l± (now = l x ) and larger / (now 
= ly) which we also present in Table [71 We conclude from this Table that: 

1) for very small I the overlap |02AS'(&)| ^ s ^ ar § e f° r a ^ 

2) and for fixed I the values of |0 2j 4s(&) I grow as l± decreases; 

3) but \02As(b)\ rapidly decreases to values consistent with zero as I — > oo, and this is so for 
any fixed b and appears to be the case for any fixed l± as well. 

We conclude that for long flux tubes on large volumes, we can safely ignore screening and 
label states as k = 2A and k = 2S. Indeed we see that it is only when I or l± are close to the 
phase transition at l c that screening is significant. Our results for k = 3A, 3M, 3S are very 
similar and the vacuum expectation value of the adjoint loop shows very similar trends. In 
practice this means that in Tables [3], [Hit is only for l/a = 16 (and l/a = 20 for some k — 3) 
that the states have needed to be extracted using the whole k = 2 or k = 3 basis (and we 
have then assigned the A, M, S labels on the basis of what component dominates the wave 
function). 

4.2 k=2A, 2S 

In the k = 2 sector we focus on the irreducible representations in / <g> /, i.e. the totally 
antisymmetric, 2A and the totally symmetric, 2S [H]. The k = 2 sector contains other 
representations, e.g. from the decomposition of / <g> / <g) / <g> /, but one expects these to have 
higher energies, and we do not consider them here. As we have remarked above, the dynamics 
appears to respect these representations very well, despite the potential mixing from gluons 
in the vacuum. Only for I ~ l c is there significant mixing. 

The lightest k = 2 flux tube is essentially pure k = 2A. We see from Table El that it 
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aE R (l = 


40; l±) 




R 


l± = 16 


l± = 32 


l± = 48 


f (k=l) 


0.2790(24) 


0.2802(28) 


0.2817(14) 


f (k=l)* 


0.527(10) 


0.522(7) 


0.507(3) 


k=2A 


0.438(4) 


0.465(7) 


0.462(5) 


k=2A* 


0.718(4) 


0.663(11) 


0.655(11) 


k=2S 


0.673(5) 


0.661(7) 


0.650(9) 


k=3A 


0.484(9) 


0.530(8) 


0.526(7) 


k=3A* 


0.774(18) 


0.719((10) 


0.732(15) 


k=3M 


0.814(7) 


0.799(24) 


0.817(13) 


k=3S 


1.14(2) 


1.116(16) 


1.154(29) 


Adj (k=0) 


0.566(6) 


0.584(10) 


0.560(9) 



Table 6: The energy, Er(1), of the lightest flux tube of length I = 40 (all I in lattice units) on 
a 40 x l± spatial volume, and with the flux belonging to the representation R. For SU(6) at 
= 171.0. 



is lighter than two fundamental flux tubes (which would also be k = 2) so this flux tube 
is absolutely stable. Its calculation therefore provides a 'benchmark' for what constitutes a 
'good' energy calculation in this paper. The energy is calculated from the correlator C(t = ant) 
of our variationally selected best trial wave-functional for the state. We can define an effective 
energy by 

ag e// (n t ) = -ln c7( ^ ) 1) (12) 

and note that if C(t) is independent of t for t > to (within errors), then this implies that it 
is given by a single exponential, C(n t )/C(0) = \c\ 2 e~ aEnt for t > to (within errors). So to 
calculate aE we need to identify a plateau in the values of aE e ff(n t ) and the quality of our 
calculation is reflected in how convincing this plateau is. 

In Fig. [T]we plot our values of aE e ff(n t ) for various values of I. We also show our final 
energy estimate in each case by the horizontal lines. We have excluded values at larger n t , 
once the errors have become larger than ~ 15 — 20% since these carry little information and 
merely clutter the plot. (In addition, at large nt the correlations within the Monte Carlo 
sequence become very long and our error estimates become increasingly unreliable.) We can 
see that we have a well-defined energy plateau for all our values of I, although the length of 
the plateau shortens as I "I since e~ Ent will disappear into the statistical noise more quickly 
with increasing t for larger E n . 

We fit these energies with the Nambu-Goto formula in eqn ffTO]) . together with a theoret- 
ically motivated 0(l/l 7 ) correction, which however plays no significant role in the fit. We 
extract the string tension a 2 a 2 A and plot in Fig [2] the values of E (l) versus I, with both 
expressed in units of the string tension. We see a very clear near-linear increase characteristic 
of linear confinement. We also see that the pure Nambu-Goto prediction appears to fit very 
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2A/2S overlap 


I l±x l t 


bl=l bl=2 bl=3 bl=4 bl=5 bl=6 


13 60 x 200 

14 60 x 200 


0.292(48) 0.344(53) 0.384(56) 0.438(58) 
0.122(25) 0.157(31) 0.187(36) 0.227(42) 


16 20 x 200 
16 40 x 200 
16 60 x 200 
16 80 x 200 
16 100 x 200 


0.172(6) 0.234(7) 0.285(8) 0.356(9) 0.491(8) 
0.053(3) 0.076(4) 0.097(5) 0.129(6) 0.206(8) 
0.036(3) 0.053(4) 0.067(5) 0.088(6) 0.136(8) 
0.034(2) 0.048(3) 0.061(4) 0.081(5) 0.122(6) 
0.032(3) 0.047(3) 0.059(4) 0.076(5) 0.116(7) 


20 16 x 200 
40 16 x 200 
60 16 x 200 
80 16 x 200 
100 16 x 200 


0.071(2) 0.125(3) 0.175(3) 0.259(4) 0.416(3) 
0.001(1) 0.002(1) 0.009(1) 0.035(2) 0.189(2) 
0.000(1) 0.001(1) 0.000(1) 0.006(1) 0.090(2) 
0.001(1) 0.001(1) 0.001(1) 0.001(1) 0.041(2) 
0.000(1) 0.001(1) 0.001(1) 0.000(1) 0.020(1) 


20 70 x 120 
24 48 x 60 
32 40 x 48 
48 48 x 48 


0.010(5) 0.015(8) 0.019(12) 0.028(16) 0.037(19) 
0.002(4) 0.001(5) 0.003(7) 0.011(9) 0.022(13) 
0.001(4) 0.003(3) 0.000(6) 0.003(7) 0.001(11) 0.107(17) 
0.001(3) 0.001(3) 0.003(5) 0.001(6) 0.004(5) 0.011(5) 



Table 7: The modulus of the normalised overlaps |02As(bl)| of blocked Polyakov loops in the 
2A and 2S representations, for blocking levels bl, as defined in eqn( TlTi) . On lattices of various 
sizes (shown in lattice units). 



well. 

Uniquely for the absolute ground state the expansion of the Nambu-Goto prediction for 
the energy Eq(1) in powers of l/l 2 a converges right through the range of I where we have cal- 
culations; indeed all the way down to ly/a = n/3 ~ 1.1 < l c ^a. This provides an opportunity 
to test not just the resummed Nambu-Goto expression, but the individual power correction 
terms predicted to be universal [3). To do this we normalise E (l) to the leading al piece, so 
that we can readily expand the scale, and compare to various 'models' for E (l). This produces 
Fig. [31 Here we see that the free string expression is good all the way down to ly/a 2 A ~ 2 
which is close to the deconfining length, l c , indicated by the vertical red line. And we note 
that a 0(l/l 7 ) correction can describe the deviations from Nambu-Goto for 1\J<Jia < 2. How- 
ever we also see that including just the leading universal correction, i.e. E (l) = al — n/6l 
[U, 12], is indistinguishable from Nambu-Goto within the errors in the range of I where the 
latter well describes E (l). However if we only include a linear al piece, then this does not fit 
at all. Thus we have a quite accurate confirmation of the presence of the universal Liischer 
correction, but not really much more than that. The reason for this is that the universal 
corrections to E (l) have small coefficients, since they represent just the zero-point energies 
of the string fluctuation modes, which indeed is why the expansion converges down to small 
I. (One can do better with the fundamental flux tube [6], since a is smaller there, and it is in 
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that case that one may realistically hope to pin down all the universal corrections.) 

It is worth quantifying how well we can constrain the Luscher correction with the k=2A 
ground state. We find 

7T 

aE 2A {l) = ~ c eff-^j ; °eff = 1-05(15) for ly/a 2 A > 2.5, (13) 

which is a usefully accurate test of this universal coefficient. 

We now turn to states with non-zero longitudinal momenta. In Fig. [2] we also plot the 
ground state energies for the lowest two non-zero momenta along the /-torus, p = 2ir/l and 
p = 4-k/I. We find that there is a unique such state for p = 2n/l and it has P = — . For 
p = 4tt/1 we find two apparently degenerate ground states, one with P = + and one with 
P = — . All this is just as expected from Nambu Goto where the p = 2n/l state has one 
phonon, and hence P = — , and the p = Air /I ground states have either one phonon carrying 
the whole momentum, with P = — , or two phonons sharing the momentum, and hence P = +. 
We also show in Fig.|2]the ground state energy of two (non-interacting) fundamental flux tubes 
of length I carrying the same total momentum. We see that this state always has a higher 
energy than that of the corresponding k = 2A flux tube showing that the latter is indeed 
stable. 

Since the only parameter in Nambu-Goto is the string tension, which is obtained by fitting 
the p = state, the Nambu-Goto predictions shown for p = 2tt/1 and p = 4tt / Z have no free 
parameters. It is therefore remarkable that the agreement is so precise and extends to our 
smallest values of I. Of course some of the energy comes from p 2 and so it is useful to perform 
a comparison with this subtracted. We therefore define the quantity: 

AE\q, I) = E 2 (q; I) - \l) - ^ * ^ 4vra(iV L + N R ), (14) 

using eqns([H]) and (ITU|) . This exposes the excitation energy predicted by Nambu-Goto. We 
plot the ratio AE 2 (q, Z)/47r<r in Fig. HI We see that the integer-valued contribution of the 
excitation energy is very accurately confirmed for all I, even for very short flux tubes which 
certainly do not 'look like' thin strings. This is something that we have already observed 
for fundamental flux tubes [B] but here we know that the flux tube is a bound state with, 
therefore, some extra internal structure. From the comparison in Fig. |5] between the k = 2 A 
energy and that of two free k = 1 flux tubes, we infer that the binding energy is not very 
large, so that at small / the k = 2A flux tube will be a 'blob' rather than a 'thin string'. It 
is therefore remarkable that its excitation spectrum should be so precisely that of a free thin 
string. 

In Fig. H]we also show what happens if one excludes the zero-point energy from the Nambu- 
Goto formula. We see a very visible shift for both p = and p = 2ir/l. (It would be pointless to 
go to higher p since the errors are too large there.) For p = this is just another presentation 
of our result in eqn fflB"]) . however it is interesting to see that the p ^ spectrum also reveals 
the presence of this zero-point energy. We do not quantify it further because it would add 
little to eqnflSP. 
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To assess the significance of these results for p ^ it is worth stepping back and asking 
what we might expect if we make no assumption at all about the relevance of stringy fluctu- 
ations. We would expect on general grounds that the absolute ground state of the flux tube 
is intrinsically translation invariant in the direction of the flux tube, so can only have p = 0. 
Thus the non-zero p has to be carried by some additional excitation. Let us suppose that this 
is some particle of mass m. Then neglecting any interaction between this particle and the flux 
tube, the energy of the combined system is 

E(l;p) = E gs (l) + (m 2 +p 2 ) 1/2 ! P=^ (15) 

where E gs (l) is the (observed) energy of the absolute ground state. To decide whether this 
model has any plausibility, we plot E(l;p) for the massless case, m = 0, in Fig.[2]as the dashed 
lines. We see that these are very close to the Nambu-goto predictions and could provide a 
good first approximation to the observed spectrum. It is therefore interesting to ask how this 
constrains the value of m. So we calculate m using eqn( ll5]) at each value of I for p = 2-k/I, 
since these p^O energies are the most accurate, and average the results for I > Iq, for various 
choices of Iq. The result, in units of the string tension, is shown in Fig. [5j (We also show 
the similar result of a similar analysis applied to flux tubes in the k = 3 A representation.) 
Roughly speaking this tells us that m 2 /<7 2a < 0.1(1). This is to be compared to the known 
value of the mass gap in the SU(6) gauge theory [231 IMj which is m^/cr 2a ~ 13. Thus this 
'particle' cannot be an excitation in the bulk space-time, and must be an excitation that lives 
on the flux tube. In that case the obvious candidate is a massless stringy mode of the kind 
described by the Nambu-Goto free string model. Note that this of course means that the 
relationship in eqn ffTol) is not the correct one. Note also that although eqn (fT5"j) is, numerically, 
very close to eqnQ for states where the massless phonons are either all right or all left movers, 
this is no longer the case when both right and left movers are present, e.g. the first excited 
p = state. As it happens, we shall shortly see that, although this state is badly described by 
the extension of eqn (fl5]) . it is also badly described by Nambu-Goto. However in the case of 
fundamental flux tubes, studied in [5J, one finds that Nambu-Goto works well for I not very 
small, and thus eqn flTSl) would be strongly disfavoured. In addition a state with a ground 
state p = flux tube and an additional particle would not couple to our operators as N — > oo 
in contrast to what one observes for the states with p ^ 0. Our purpose in considering this 
simple model was to establish, in a pedestrian way, that one must look to massless modes 
living on the flux tube for the origin of the observed spectrum. 

We turn now to the spectrum of excited states with p = 0. We plot, in Fig [HI the four 
lightest P = + states, and the two lightest P = — ones, as well the predictions of Nambu- 
Goto for the lowest few energy levels. (We also plot some higher excitations for / = 32a and 
I = 52a, which we shall return to shortly.) In Nambu-Goto the ground state, with no phonons, 
is non-degenerate, with P = +, as is the first excited energy level which has one left and one 
right moving phonon with momenta p = ±2ir/l. The next energy level has four degenerate 
states with the left and right moving phonons sharing twice the minimum momentum. Since 
this can be carried by one or two phonons, two of these states have P = + and two have 
P = — . If the 2 A flux tube states were close to Nambu-Goto, as they turn out to be for 
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the case of fundamental flux, we should find our calculated energies clustering closely about 
the lowest three Nambu-Goto energy levels. While we do indeed observe in Fig [6] that the 
lightest two states do have parity P = +, and the next two P = + states are roughly in the 
same energy range as the lightest two P = — states, we see nothing like the (near) degeneracy 
predicted by Nambu-Goto. There is some evidence that the first excited P = + state and 
the lightest P = — state approach the appropriate Nambu-Goto levels, and that the second 
lightest P = — state agrees with the Nambu-Goto prediction for all but the smallest values 
of /. However the observed excited states are , in general, far from showing the Nambu-Goto 
degeneracies and are far from the Nambu-Goto predicted energies, even for the largest values 
of ly/a. While the first excited state appears to clearly approach the string prediction, even 
here it would be useful to have some further evidence that it is asymptoting to that curve and 
not just crossing it. It is useful to recall that for the fundamental flux tube [6], the convergence 
to Nambu-Goto was rapid and unambiguous (albeit not as rapid as for the ground state). The 
messiness of the picture in Figgis of course what one would have naively expected for such a 
bound state flux tube, and the real surprise is the precise stringy behaviour we have observed 
for the lightest states with non-zero momenta. One significant difference with the latter is that 
here the states are generally well above the threshold for decay. The lightest asymptotic decay 
products will be two fundamental flux tubes with equal and opposite transverse momentum. 
The energy of the threshold, corresponding to zero relative momentum, is shown in Fig |6] and 
one is tempted to note that the deviation of the first excited state from Nambu-Goto decreases 
as the phase space for decay decreases. We also show the energy of a decay state composed 
of a glueball and a ground state k = 2A flux tube. We see that this is quite high and, in 
addition, such a decay will be large- N suppressed. 

To provide some more context for these states, we have also shown in Fig [6] the next 
6 P = + and 5 P = — states for I = 32a and I = 52a (slightly shifted in I for clarity). 
The number of states is motivated by the fact that the next Nambu-Goto energy level has 5 
P = + and 4 P = — degenerate states, so we are also including at least one state, for each P, 
that will approach a yet higher energy level as I — > oo. (But note that the extraction of the 
energies can be ambiguous for these massive states.) The main message, considering all the 
I = 32a, 52a states, is that there is no visible clustering in the energy of the states that might 
suggest that they are converging to the Nambu-Goto energy levels, except for the absolute 
ground state and perhaps the first excited state, both of which are P = + at our largest 
value of I where a clear gap has opened between them and the P = — states - as expected in 
Nambu-Goto. For the first excited state there is a residual ambiguity: is it the first excited 
state at lower I that asymptotes to the Nambu-Goto level as I — > oo, or is it perhaps the 
second, with the first 'crossing' that level somewhere between I = 52a and I = 64a ? In fact 
our analysis in Section [5] will address and resolve this issue. What we see in Fig[6j particularly 
for I = 52a, is very much a continuous distribution of excited states without any obvious 
level structure. This makes it hard, for example, to know whether the near-coincidence of the 
second P = — energy with the Nambu-Goto prediction is in fact significant, or merely the 
chance result of this near-continuous distribution of states. What is clear from Fig O is that 
we are very far from the values of I where the Nambu-Goto spectrum might become a good 
first approximation for these states. 
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The fact that these excited p = states are generally well above their decay thresholds 
raises some questions. The most important is how confident can we be that we have extracted 
their 'energies'? If the decay width is very small, the propagator should have a pole in the 
complex energy plane very close to the real axis, and we would expect correlators designed for 
stable states to behave just as they do for a stable state, within the finite statistical errors, i.e. 
we should see an effective energy plateau that is lost in the statistical errors at larger t before 
deviations from the plateau become visible. We show the effective energy plateaux for the 
first excited state in Fig. [71 We see that for large / these plateaux are unambiguous and not so 
different from those of our stable ground state in Fig. [U As I decreases, however, the apparent 
plateau shifts to larger t and becomes increasingly ambiguous. This is very different to what 
we observe in Fig. [TJ The likely reason for this is that the phase space for the decay of the first 
excited state grows as I decreases (as we scan infer from Fig [6]), and so presumably does the 
decay rate. So it is interesting to perform a different analysis, at the smallest values of I, that 
attempts to take this finite decay width into account. The relevant asymptotic states in this 
energy range are those composed of two (unexcited) fundamental flux tubes with equal and 
opposite transverse momenta. (Flux tubes with longitudinal momenta have larger energies.) 
Obviously if we performed a variational calculation with a complete basis of k = 2 operators, 
then these are the states we would obtain. However the operators we actually use are all of 
the form Tr 2 A{^ P } oc Tr/{/ p } 2 — Tr/{/p} with l p some winding operator. The Tr/{/ p } 2 piece 
represents two fundamental flux tubes at zero spatial separation, which can be re-expressed 
as a sum over all relative momenta. However the projection onto any such state with given 
momentum will be very small, so a variational calculation performed at t — a, as ours is, will 
not pick out these states. However the overall projection onto all these states is not small, and 
a heuristic procedure is to perform a fit to the correlation function that is in terms of these 
asymptotic scattering states, but with an amplitude that encodes a slightly unstable state. 
We choose, again heuristically, a Breit-Wigner form. So we fit to: 

C(t) = J2\cBw\ 2 e- Eip)t ; E(p) = 2E f (p), \c BW \> = ° (16) 

g^J (E - E ) + (E L ) 

where Ef is the lightest energy of a fundamental flux tube with transverse momentum p, c is a 
constant fixed by normalisation, Eq is the real part of the pole energy and To the (full) width. 
We either sum over a discretisation of the momentum integral, or use the transverse momenta 
dictated by the size of the transverse torus. (In practice it does not matter which we use.) In 
Fig M we display the values of aE e /f for I = 16a, on a blown- up scale, and display different 
fits. The red line arises from a conventional fit with an excited state in addition to the desired 
lightest state. Here the lightest state is at aE = 0.25, and the heavier one is at aE = 0.49, with 
relative probabilities 15% and 85% respectively. With such a low overlap, we can have little 
confidence in the robustness of this lightest state. The alternative fit based on eqn( fl6|) is shown 
by the solid black line and corresponds to aE = 0.475 and r = 0.065. (The dotted black line 
corresponds to a sum over scattering states with uniform probability.) We see that the value 
of the energy is close to but larger than E e ff(t = a) whereas a search for a large- 1 plateau 
leads (as in our first fit above) to a much lower result. This is characteristic of such fits. We 
note that it is no coincidence that in our first, conventional, fit the dominating 'excited' state 



17 



at aE = 0.49 is close to our Breit-Wigner pole in the alternative fit based on eqn( [16]) . This 
gives us confidence that this is most likely the actual energy of this unstable excited state. We 
note that applying such a procedure would raise the energy estimate significantly closer to the 
Nambu-Goto prediction. For example, a similar analysis at I = 20a would give aE p = 0.485, 
with T p = 0.050, rather than the value aE ~ 0.38 from a plateau estimate, and this would 
approximately halve the discrepancy with Nambu-Goto. The effect is even more marked at 
I = 16a. Clearly what we need is sufficient statistical accuracy to distinguish between the two 
different n t — > oo values of aE e ff(n t ) in Fig[HJ Moreover it would be useful to see the stability 
of such an analysis to the presence of a second heavier excited state (which surely contributes 
at some level). Nonetheless, while we cannot be definitive on this, it is plausible that where 
the apparent plateau is indistinct because it is at large n t , and in addition the state has a 
large phase space to decay, the actual energy of the 'resonant' flux tube is much closer to the 
value of aE e ff(t — > 0) than to aE e ff(t — > oo). In the present case this would suggest values 
for the first excited state that are closer to the Nambu-Goto prediction at small I than our 
conventional estimates shown in Fig So it is not possible for us to be certain how much of 
the large apparent deviation from Nambu-Goto is due to the extra modes associated with the 
internal structure of the k = 2A flux tube, and how much is a consequence of the fact that 
these states are unstable. 

Two remarks. The first is that none of the above caveats apply to the ground states with 
p 7^ shown in Fig. [2j Here the effective energy plateaux (which we do not show) start at 
small t and typically become increasingly well-defined as I decreases. The second remark is 
that one might wonder if some of the apparent downward drift in E e ff(n t ) that we see at large 
n t in Fig. [7] is not due some small admixture of the ground state in our variationally estimated 
excited state wave function. Since our variational ground state wave function has a typical 
overlap onto the ground state of ~ 0.985(15) (which can be inferred from the E e ff(n t ) values 
shown in Fig. [T]) we can estimate the maximum such contribution to the excited E e ff(n t ) in 
Fig. [7J and it turns out to be invisible for l/a > 24 (at our level of accuracy) and only possibly 
becomes visible for nt > 15 for l/a = 16, 20. That is to say, it is essentially irrelevant here. 

We turn now to flux tube states obtained by performing calculations with operators pro- 
jected onto the k = 2S representation, i.e Tr 2 5{/ P } oc Try{Z p } 2 + Tr f{l^} with l p some winding 
operator. (For the / = 16a p = ground state we obtain a cleaner variational state by using 
the full 2A © 25" basis, and that is what we show here. The admixture of 2A is small and so 
it still makes sense to label the state as 2S, as we do.) We know that these will be heavier 
than the corresponding k = 2A states [161 E] and so we expect all of them to be unstable 
as well as having larger statistical errors. In Fig. ED we show the ground states with the low- 
est longitudinal momenta. We also show the energies of the lightest decay products in each 
case. Just as for k = 2A the energies are remarkably close to the Nambu-Goto predictions, 
as emphasised by comparing the actual excitation energies in Fig. [TUJ In Fig. [11] we show 
the effective energies for the absolute p = ground state. We also indicate the expected 
decay thresholds on the right side of the figure. It seems clear that E e ff(n t ) does possess 
extended plateaux very different from the decay thresholds in the lowest / cases where we 
have accurate results to large t. So while the quality of the calculations is markedly inferior 
to the k = 2A case, we have confidence in our extraction of the energies plotted in Fig. [9j The 
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situation with the excited p = states is however much worse and we are unable to extract 
the corresponding energies. Our problem is illustrated by Fig. [12] where we plot the effective 
energies for the 'state' selected by our variational procedure as the first excited p = state. 
We cannot identify a plausible energy plateau for any value of /, and E e ff(n t ) is consistent 
with a decrease towards the decay thresholds shown. In Fig. [13] we repeat the exercise in 
Fig [8] now for the l/a = 16 k = 2S flux tube. The fit using eqn fll6l) works very well, and 
corresponds to an energy E = 0.58 and a width T = 0.1. The two exponential fit is less 
convincing and corresponds to energies 0.205 and 0.595 with overlaps squared of 0.1 and 0.9 
respectively. This begins to point rather unambiguously to an energy estimate of E ~ 0.58 
and hence Ej^/o2s ~ 4.6 at l^/o2s ~ 2.1. We note that this is below, but not far below, the 
Nambu-Goto prediction. A similar conclusion follows for l/a = 20. It is thus plausible that 
this unstable first excited state is indeed quite close to Nambu-Goto although this would be 
far from apparent using a conventional analysis. 

4.3 k=3A, 3M, 3S 

In the k = 3 sector we focus on the irreducible representations in / <g> / <g) /, which are the 
totally antisymmetric, 3A, the mixed, 3M, and the totally symmetric, 3S. We know from 
earlier work [TBI E] that the corresponding string tensions are very close to the predictions of 
Casimir scaling (see also Table EJ) and so, as we shall see, the ground 3^4 states are stable, the 
3M states nearly so, and the 35* states are highly unstable. 

In Fig. [14] we plot the lightest energies of k = 3A flux tubes with longitudinal momenta 
p = 0, 2tt/1, 4n/l. Just as for the corresponding k = 2 flux tubes, we see excellent agreement 
with Nambu-Goto all the way down to I ~ l c . The relevant asymptotic decay states are 
not just 3 fundamental flux tubes, but also a stable k = 2 A flux tube with a fundamental 
one. The latter is lighter and the thresholds for both are plotted as the black lines in Fig. [TT1 
demonstrating the stability of the k = 3A states. As we see in Fig. [15] for the absolute ground 
state, we have very well defined energy plateaux, again just as for the k = 2A flux tubes. 

As for the k = 2A case, it is worth quantifying how well we can constrain the Liischer 
correction with the k = 3A ground state. Here we find 

7T 

aE 3A (l) = a 3A l - c eff — ; c e// = 1.11(11) for l^/a 3A > 2.3, (17) 

which is again a usefully accurate test of this universal coefficient. 

We turn now to the lightest excited states in the p = sector, as displayed in Fig. [TB] 
Comparing to Fig. IB] we see that the phase space for the first excited flux tube to decay is 
smaller here and indeed at larger I it is stable. This is perhaps why its energy, particularly at 
small I, is closer to Nambu-Goto than in the k = 2 A case. And also why the effective energies 
displayed in Fig. [T71 show clear plateaux even for l/a — 16, in contrast to the k = 2 A case in 
Fig. [7] (Note that for l/a = 16 we use our full k = 3 basis, which means that state includes 
a very slight admixture of k = 3M and k = 35.) The decay thresholds are indicated on the 
right hand axis of Fig. [T71 and it is clear that the low-Z plateaux take very different values. 
This provides us with a quite clean example of an excitation of a bound state flux tube where 
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we can ignore the (slight) instability of the state. It is therefore interesting to compare this 
to the corresponding excitation of the fundamental flux tube in Fig. 19 of |Bj. We see that 
the deviation from Nambu-Goto is indeed very much larger here, and this must be due to the 
bound state structure of this flux tube. We note that a similar analysis applied to the first 
P = — excitation with P = 2n/l leads to very similar conclusions. 

We turn now to the heavier k = 3M states. We plot in Fig. [18] the ground states with 
the lowest longitudinal momenta. Once again these particular states agree very well with 
the Nambu-Goto predictions. However we see that they are now slightly above the decay 
threshold and so will be unstable but apparently not enough to affect the extraction of, for 
example, the absolute ground state as we see in Fig. [191 (Again we use the full k = 3 basis 
for the 1 1 a = 16 ground state.) However the p = excited states are very unstable and we 
are unable to identify useful plateaux. 

The k = 3S states are much heavier and we can only estimate energies for the ground state 
p = 0, 2ix/l states, as shown in Fig. [2U1 Again we see rough agreement with Nambu-Goto, but 
now the decay phase space is large - becoming very large for large I. We show the effective 
energies for the absolute ground state in Fig. [2TJ While the plateaux at lower I are quite clear 
and are far from the decay thresholds (indicated on the right hand axis), this is not the case 
at the largest values of /. (Indeed we do not even attempt to extract an energy for I /a = 48.) 
In the latter cases, while the motivation for our energy estimates should be apparent, it is not 
necessarily convincing. Nonetheless the usual agreement with Nambu-Goto for such states at 
smaller I is remarkable. 

4.4 adjoint 

The adjoint flux tube appears in /<g) / and should couple to operators Tr a djl p = |Tr// p | 2 — 1, if 
indeed it exists. There is some evidence from the calculation of adjoint potentials that it does 
indeed exist and that the adjoint string tension satisfies approximate Casimir scaling (see e.g. 
[j~9] and references therein). Such a flux tube can be screened down to the vacuum by gluons, 
but this is suppressed by 1/N 2 , and is in fact negligible except for finite volume effects. The 
latter can either arise if I is small, i.e. I ~ l c , or if we consider blocked/smear l v operators 
that extend around the transverse torus. In practice we always include such highly smeared 
operators in our calculations, since they (slightly) improve the overlap onto the ground state 
of the adjoint flux tube, and we therefore explicitly subtract vacuum expectation values in 
our correlators. 

An adjoint flux tube whose string tension satisfies approximate Casimir scaling will in 
general be heavier than a pair of fundamental anti-fundamental flux tubes and can therefore 
decay into these. (Here there is no large- N suppression.) Just as with unstable fc-strings, 
the important question is whether the adjoint flux tube is nearly stable, so that conventional 
methods for extracting the energy can be used, or not. We shall be careful to establish whether 
this is so or not. 

In Fig. |22]we plot the energies of the lightest adjoint flux tubes with longitudinal momenta 
p = 0, 2tt/1, Ait /I. As usual the p = Nambu-Goto fit fixes the string tension a 2 a a dj, and then 
the Nambu-Goto predictions for p ^ are parameter- free. We observe that, again as usual, 
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these predictions are remarkably well satisfied all the way down to / ~ l c . The decay thresholds 
are indicated and we see that the decay phase space is small, raising the hope that the decay 
widths will be negligibly small. Of course the statistical errors are quite large here so it is 
worth extracting the 'excitation energy' as defined in eqn (H4p to see how well that is being 
determined. As we see from Fig. [23] the modes carrying momentum are indeed unambiguously 
the wave-like modes of a thin relativistic string. 

In Fig. [24] we plot the effective energies for the absolute ground state. (Energies shifted for 
clarity.) Horizontal red solid lines indicate our plateaux estimates, including errors. For small 
and medium I these are well determined, but for the largest values of I the states are very 
massive and we quickly lose the signal as we go to larger n t . Hence the generous error estimates 
in these cases. For comparison we plot the // threshold energies as horizontal dashed lines 
(also as points on the right hand axis). These are quite close to the plateaux, especially at 
small I. So we blow up the scale for the latter states in Fig. [25] A characteristic feature of 
effective energies is that once the error gets large, the estimate of that error becomes unreliable. 
This applies to the large n t decrease or increase in aE e ff(n t ) that we see in Fig. [25] Since 
our correlators are diagonal, an increase would violate positivity, and so must be statistical. 
There is therefore no reason to take the decreases any more seriously. Given these remarks, 
we can see that the l/a = 16 plateau estimate is consistent with the decay threshold, while 
for l/a > 20 (and unambiguously for l/a > 20) the plateaux is well above the threshold. We 
conclude that the adjoint flux tube does indeed exist as a distinct and nearly stable 'bound 
state'. 

On the other hand we cannot identify well-defined excited states with p = 0. These would 
have a very large phase space for decay into // flux tubes, so this is not unexpected. They 
are presumably analogous to broad resonances, and will be equally difficult to identify. 

4.5 84 and 120 

In the / <g) / (g) / sector of SU(6), the irreducible representations with the smallest Casimirs 
and, we can assume, the smallest string tensions, are the 84 and 120 . (See the Appendix.) 
Here we shall study flux tubes carrying flux in these two representations. 

Such flux tubes can mix with single fundamental flux tubes, but this is large- N suppressed 
and given our experience with the adjoint flux tube, we shall (usually) ignore this possibility. 
However the decay/mixing with 3 (ant i) fundamental flux tubes is not large- iV suppressed. And 
neither is that with a k = 2A and an antifundamental, which is even lighter. In Fig. [26]we plot 
the energies of the ground state 84 flux tubes with longitudinal momenta p = 0,2^/1,4^/1. 
The Nambu-Goto predictions are shown as solid red curves, with the decay 3/ and 2A + / 
thresholds indicated by the black curves. As usual we extract the string tension from the 
p = fit so that the p ^ predictions are parameter free. We observe that the agreement 
is, once again, remarkably good for p = 2n/l and quite good for p = 4ir/l, where however 
the states are very massive and it becomes difficult to identify plausible plateaux. The string 
tension is comparable to that for k = 3M (see Table [2]) as is the phase space for decays. So it 
is no surprise that, just as for 3M, we are unable to obtain energy estimates for p = excited 
states. The ground state however has reasonably clear energy plateaux, as we see in Fig. [27] 
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at least for 20 < l/a < 40. For l/a > 44 the effective energies are large and disappear rapidly 
into the statistical noise as n t increases, making plateau identification increasingly subjective. 
For l/a = 16 we see no plateau, and here we see that E e ff(n t ) decreases well below the decay 
thresholds shown and appears to be asymptoting to a large- N suppressed single / admixture. 
That this should only occur for our shortest flux tube, l/a = 16, is consistent with our earlier 
observations about the finite volume effects displayed in Table [7J 

In Fig. [28] we plot the ground state energies of flux tubes in the 120 representation for 
p = 0, 27r/Z. The 120 string tension, which we obtain by fitting the p = values, is almost as 
large as the k = 3S one, and so it is no surprise that just as in that case we have no useful 
results for p = Air /I or for any excited states. Indeed even the p = effective energy plateaux 
are difficult and ambiguous to identify in this case. 

Finally we remark that we have also performed some matching calculations in SU(3) at 
(3 = 40.0, which corresponds to about the same lattice spacing. The corresponding / <g> f <g> f 
irreducible representations are the 6 and 15. In both cases the energy plateaux are more 
ambiguous, particularly where we compare the 15 with the 120 of SU(6). This may be due to 
the fact that certain mixings and decays are less suppressed for SU(3) than for SU(6). 

5 Excited states: massive or stringy? 

One of our motivations for studying flux tubes in higher representations is that we expect 
such bound states of fundamental flux tubes to have a low-lying excitation spectrum that 
contains clear signatures of the binding scale. This should provide an interesting contrast to 
the low-lying spectrum of fundamental flux tubes which, unexpectedly, shows no sign of the 
excitation of the massive modes that one would expect to be associated with an 'intrinsic 
width' for the flux tube. While one might question the existence of such an intrinsic width, 
the existence of a non-zero binding in the case of, say, the 2A flux tube is unambiguous. This 
would, most simply, reveal itself in extra excited states, representing massive rather than the 
usual stringy massless modes. Our cleanest spectra in this paper are for k = 2A and k = 3A 
so we shall focus on these. So does the k = 2A p = spectrum shown in Fig. El reveal any 
massive modes that are additional to the stringy excitations which, at large /, tend to the 
Nambu-Goto curves? (The same observations apply to the k = 3A spectrum.) Since the low- 
lying excitation spectrum of fundamental flux tubes appears to contain only stringy states 
and no massive modes, it is interesting to compare our k = 2A spectrum to the fundamental 
one shown in Fig. 12 of [6]. The immediate question this comparison raises, as pointed out in 
our earlier study of k = 2 flux tubes in [7], is whether the first excited k = 2 A state might 
be a massive mode, with the second excited state being the first excited stringy mode and 
the next two P = + excited states eventually tending to the second Nambu-Goto level? (We 
have not shown higher P = + excited states in Fig. [61 but they are there.) Or it might be 
that the large deviations from Nambu-Goto are largely driven by the 'unstable' character of 
these flux tubes, and that otherwise the first excited state is much like the fundamental one. 
However this possibility appears to be contradicted by our results in this paper for the much 
more stable k = 3A states, plotted in Fig. [161 which show similarly large deviations from the 
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Nambu-Goto predictions. Or again, it might be that we are seeing here the mixing of modes, 
enhanced by the existence of intermediate states that are not far from threshold. This could 
be the mixing of nearby stringy modes, or of a stringy mode with a massive mode - which 
would also imply the presence of an extra mode. 

So we want to ask if the first excited states in the k = 2A and k = 1 cases are the 
'same' or not. It is of course not possible to answer this question unambiguously, and we 
choose to address it in the same way as we did in [Tj. The idea is that if this state is 
indeed an approximate Nambu-Goto-like string excitation then we would expect its wave- 
functional to have the appropriate 'shape'. What that 'shape' should be, in terms of our 
highly blocked/smeared link matrices, is not at all evident, but it is something we do not need 
to know because we can simply compare it to the wavefunctional of the first excited k = 1 
state, which we have good reason to think of as being stringy. 

The way we make this comparison is as follows. Let {<pi]i = 1, n„} be our set of 
winding flux tube operators, with P = + and p = 0. These operators are group elements, 
not yet traced, and may be in any representation of SU(iV). Suppose the flux is in the 
representation R. When we perform our variational calculation over this basis, we obtain a 
set of wavefunctionals, $> R , which are an approximation to the corresponding eigenfunctions 
of the Hamiltonian. Unfortunately we cannot simply compare R = f and R = 2A states by 
calculating their overlap: it will vanish because of the center symmetry. So instead we proceed 
as follows [7] . We write the wavefunctionals as linear combinations of our basis operators: 

$ R = E bliCR^R^) = J2 ^:T4(0>) (18) 

i i 

choosing the coefficients cr^ to satisfy the normalisation condition 

(TrJ(0 i (O))Tr^(0 i (O))) = 1 (19) 

so as to ensure that a comparison of the coefficients b Ri between different representations R 
can be meaningful. The idea is that the coefficients h Ri encode the 'shape' of the state corre- 
sponding to the wavefunctional, because they multiply the same operators, albeit in different 
representations, and with a common normalisation. So making the simple substitution 

$Li = E^r' 2A (0 4 ) — ► = XX/BW ( 2 °) 

i i 

we can compare our excited k = 1 and k = 2A wavefunctionals by comparing with the 
fundamental wavefunctionals, This we can do by calculating the overlap 

O = ^/^a) (21) 

<n ($"'t / $/) 1 /2(|»«t| ) n A) l/ 2 

(all operators at t — 0) which we assume provides us with a measure of the similarity between 
the original state an d the state $j . 
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Even if one accepts this method of comparison, there are some important caveats. The 
variational calculation is performed over a limited basis, so the $^ are only approximate 
energy eigenfunctionals. And the level of approximation will generally be different for different 
representations (and states). Thus the comparison is inevitably approximate. Again, we note 
that the operator basis (pi varies with I (in lattice units). So we perform the comparison of 
/ and 2 A states at the same I. However ideally we should also compare at the same string 
tension i.e. at different lattice spacings such that ay/a 2 A = a 'y/ a f an d hence different f3. 
Because of the additional costs we have not done so here, and this also makes the comparison 
approximate. 

Given the approximate and heuristic nature of this method, we need to test it in a case 
where we are confident that we know the answer. This is the case for the absolute ground 
state. So in Fig. [29] we display the above overlaps of the variational ground states of the 
k = 2A and adjoint flux tubes onto the lightest 20 fundamental variational eigenfunctionals, 
all on l/a = 32 lattices. In Fig. [301 we do the same on a l/a = 64 lattice. The result is 
clear-cut, both for stable and unstable flux tubes, and for both lengths: we observe that the 
method works very well in producing an almost exclusive overlap onto the / ground state. 
This is in fact representative of all our results for the absolute ground state, even where the 
state is unstable, and this gives us some confidence in this method. 

We turn now to the k = 2A p = first excited state. In Fig. [31] we show the overlap 
of onto the lowest for l/a = 32 and l/a = 64 lattices. In Fig. [32] we do the same 
for the corresponding k = 3A state. While the largest overlap is indeed on the first excited 
fundamental state, there is also a small but visible overlap onto the 2nd excited stringy / 
state, $j =2 , and this is very similar for the k = 2 A and k = 3 A states and for both flux tube 
lengths. While the comparison is not as unambiguous as for the ground state, it is hard to 
avoid the conclusion that this state is definitely not some new massive mode excitation. Rather 
it appears to be largely the first excited stringy mode, with a modest admixture of the second. 
The shift in energy away from Nambu-Goto might be largely the result of this mixing. We also 
note that the mixing appears to become smaller as / increases from l/a = 32 to l/a = 64 and 
the energy approaches that of the Nambu-Goto prediction. We remark that all this confirms 
that the first excited state at lower I does indeed asymptote to the first excited Nambu-Goto 
energy level, and does not cross the latter somewhere between l/a = 52 and l/a = 64 - an 
alternative possibility that we discussed earlier, in Section I4.2[ when considering Fig. [6] 

It is interesting to contrast this with what one finds for the 2nd and 3rd excited states in 
the k = 2A, 3A and p = 0, P = + sectors. Typical examples are shown in Figsj33] and [3H 
Here it is hard to draw any conclusion. While the dominant overlap is onto the corresponding 
fundamental excited state, there is a large projection on other states as well. It certainly 
appears possible that some new massive mode either dominates or is mixed into one or both 
of these states. 
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6 Discussion and conclusions 



In this paper we have calculated the low-lying spectrum of closed flux tubes in various repre- 
sentations, with the length of the flux tube stabilised by closing it around a spatial torus. We 
had several motivations for this study. 

One is to compare the resulting spectrum to simple effective string actions, just as we 
did in our earlier work on fundamental flux tubes [6]. Since higher representation flux tubes 
can be thought of as bound states of (ant i) fundamental flux tubes, the massive excitation 
modes associated with that binding should leave a signature in the spectrum. In the case of 
fundamental flux tubes we found no trace at all of massive non-stringy modes and our hope 
was to find something different here. 

Of course only a few of these flux tubes are stable against decay and are real bound states. 
Recall that only some decays are large-iV suppressed by large-iV counting arguments. For 
example the decay of an adjoint flux tube to the vacuum (plus glueballs) is suppressed at 
iV = oo, but its decay into a pair of noninteracting fundamental and anti-fundamental flux 
tubes need not be, since Tr ad j/ p = Trf/ p Tr f /^ — 1. (Although the dynamics may of course 
suppress such decays.) Most higher representation 'flux tubes' are unstable at large N, in this 
sense, and it is not a priori clear if they exist in the same way as unstable 'resonant' particle 
states exist. While there is evidence for the existence of such flux tubes when attached to 
appropriate sources [TTJ] , their stability in that case is usually ensured (for the relevant length 
scales) by the fact that screening the sources by gluons costs extra energy. Our closed flux 
tubes are not protected from being screened, and decaying, in this way. So one of the things 
we wished to learn is which of the unstable flux tubes were stable enough that one could 
analyse them by conventional methods, and which needed new methods and what those new 
methods might be. 

A closely related question is whether it makes sense to classify flux tubes according to 
the irreducible representations of SU(iV), given that the vacuum contains adjoint gluons that 
can screen and mix the representations of sources and flux tubes. Earlier calculations have 
provided evidence that this is indeed the case for fc-strings [HI E] , and in Table [7] we have 
provided similar evidence for all the representations being considered here. Apart from finite 
volume effects (both in the transverse size and in I) the screening appears to be (almost?) 
exact. Why this should be so, at what is not a very large value of N, is an interesting puzzle 
that is being investigated by us, more systematically, elsewhere [20J. 

Our first conclusion, discussed in Section HI is that the absolute ground state and the 
lightest states with non-zero longitudinal momenta are accurately described by the free string 
expression in eqn([n]) all the way down to our lowest values of I, which are very close to the 
minimal possible flux tube length at / = 1/T C . This is clearest for the very stable flux tubes in 
the k = 2 A and k = 3A representations, which are our most accurate calculations (see Figs [2] 
and [HD but it is also the case, within larger errors, for all representations, including flux tubes 
that could be very unstable. (See Figs© Q2J E3 |22J [23 M for the 2S, 3M, 3S, adjoint, 84, 
and 120 representations respectively.) This is of course just what has been observed for the 
corresponding fundamental flux tube states in [6]. 

In the case of the p = ground state the stringy corrections are small all the way down 
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to lyfo ~ 1 (basically because they arise from the zero-point energies) and so an expansion 
of the energy in powers of l/ly/cr is convergent over our whole range of I. So it is fair to 
claim that the close agreement we are seeing with Nambu-Goto down to small values oily/ a 
is a prediction of the known universal corrections to the linear al piece [31 [121 El E3] - at 
least when the flux tube is stable. In fact, as discussed in Section I4.2[ what we are able to 
confirm, within our statistical errors, is the presence of the linear piece and, in some cases, 
the — 7r/6Z universal Liischer correction, but not really any more than that. (See e.g. Fig. [3j) 
Our most accurate spectra, for the k = 2A and k = 3A representations, allowed us to confirm 
the universal value of this coefficient at the ±10% level, which is a usefully accurate result for 
these bound-state flux tubes. 

The p ground states are another matter. Here we can confirm, in some cases very 
accurately, the excitation energy = \p\ of the massless excitation that carries the momentum 
on the background flux tube, as shown in FigsJH EH [23] for the 2A, 2S and adjoint flux 
tubes respectively. In fact our calculations are accurate enough to confirm the presence of the 
additional zero-point energy, as discussed in Section 14.21 and displayed in FigHJ Indeed we 
saw that in any non-stringy attempt to describe these spectra, the particle excitation carrying 
the non-zero momentum will have a mass that is constrained by our calculated spectrum to 
be very much smaller than the known mass gap of the bulk space-time theory [231 121] • Thus 
such a (presumably massless) excitation must exist on the flux tube rather than in the bulk, 
and will thus arise from an effective string action. 

That we observe a (near) free-string behaviour, even when the flux tube is very short, 
is more surprising for the p ^ ground states than for the p = ground state. This is 
because the expansion in powers of \jl\Ja diverges at quite large I for p ^ 0, so we cannot 
use universality arguments to predict the p / spectrum at small I in the way we could 
for p = 0. All this parallels what has previously been seen for the fundamental flux tube 
[5J. There we suggested [BJ that what these ground states have in common is that they all 
have a phonon content that is either all left moving or all right moving, so that all phonon 
subenergies are at threshold and it is plausible that the interaction of these Goldstone bosons 
will vanish there, removing at least one possible source of corrections to the free- string result. 
We note that a very recent and much more complete analysis of phonon scattering in a finite 
volume comes to a similar conclusion [TJ]. Our results here, with higher representation flux 
tubes, are consistent with this picture. We also remark that, as shown in Fig 12, a 'minimalist' 
model where the ground state with p ^ consists of a p = flux tube together with massless 
noninteracting particles sharing p, gives predictions very close to Nambu-Goto. (In contrast 
to a large discrepancy for other states that involve both right and left movers.) It may thus 
be that, to a first approximation, this part of the spectrum is independent of the model used 
(within limits). But this is only a speculation. 

By contrast other states, for example the p = excited k = 2A and k = 3A states shown 
in Figs [6J and [161 show very large deviations from Nambu-Goto, making it hard to say to what 
extent that model provides any kind of first approximation to this part of the spectrum. We 
recall that in the case of fundamental flux tubes the corresponding corrections are significant 
but small, with a rapid approach to the free-string spectrum as / 1 0- There are two obvious 
differences between the fundamental flux tube and the ones here. Firstly, here we have a 
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binding dynamics which may perturb the spectrum through its excitations. Secondly most of 
the states are not only unstable at finite N but remain so in the N — > oo limit. If mixing is 
the dominant effect, then the states will contain non-stringy components, while if instability 
is the dominant effect then they will be stringy but resonance-like. In this context we note 
that at the larger values of / in FigJT6l the lightest p = excited states become stable, in 
contrast to the case of k = 2A shown in Figj6j Nonetheless the deviations from Nambu-Goto 
are not significantly different in the two cases. This suggests that instability is not the main 
reason for the large deviations we see for the second excited state at these values of I. 

The striking difference between the simple stringy behaviour of the p ^ ground states, 
even when these are unstable, and the messy behaviour of the unstable excited states may be 
due to the fact that in the former case, unlike the latter, the phonons will not catalyse the 
decay of the flux tube because they have zero subenergies (3|. Although this may not be the 
only source of flux tube decay, it may be a significant factor. 

To make some progress in identifying the nature of the p = excited states we intro- 
duced in eqns( 1THH2"T]) a heuristic measure [7] for comparing states in different representations. 
Applied to the p = ground states it confirmed unequivocally that these are just the same 
as the unexcited flux tube in the fundamental ground state, and that this is so for all the 
representations we consider here. See for example Figs [29] and [30] for the k = 2 A and adjoint 
cases. This is as expected, and motivates the use of the measure for the more controversial 
excited states. For the first excited k = 2A p = state, as shown in Fig. [3U the result is less 
unequivocal but points to it being quite similar to the corresponding fundamental state and 
becoming more so as I f. We note that what we show in the figures is the overlap-squared, 
whereas it is possible that the energy shift contains pieces proportional to the overlap, which 
is larger and could produce a significant shift in the energy. For the second and third excited 
p = states, analysed in Figs [33 and [3U we see states that are very unlike the corresponding 
fundamental ones. This suggests that while the ground and first excited p = states in Fig. 
are indeed (mostly) stringy, the higher states may well include a large admixture with massive 
modes. 

Because nearly all our states are unstable, it is important to approach our energy estimates 
as critically as possible, and we have attempted to do so. (At the risk of being tedious.) One 
usually obtains an energy from a correlation function by calculating the effective energy, 
aE e ff(t) = — lnC(t + a)/C(t), and identifying a 'plateau' for t > t where t has to be small 
if the result is to be usefully accurate. The stable k = 2A, 3A p = ground states provide 
our benchmark for an unambiguous calculation, as shown in Figs [1] [15] As I f the energy 
increases oc / + 0(1//) and the signal, ~ exp{— aE(l)n t }, drops into the statistical noise at 
smaller t = an t , so decreasing the useful extent of the plateau. Apart from this there is no 
ambiguity in extracting energies. The p ^ ground states are similar, although the energies 
are larger so, once again, reducing the useful energy plateau. States that are unstable but 
with a small decay width should, by continuity, produce a plateau at intermediate t, which 
eventually sinks to the decay threshold. (As long as we remain within the same limited class 
of operators that are designed to project onto single flux tubes.) If the decay width is large, 
we will lose all sign of a plateau. Examples of the former are the adjoint ground state in 
Figs EH and [25] and, less stable, the k = 2S ground state in Fig. [11] A totally unstable state 
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is the first excited k = 2S p = state shown in Fig. [12j For obvious reasons we have not 
attempted to extract an energy for this state. This is to be contrasted with the first excited 
k = 2A p = state which we show in Fig. [71 While there is no plateau for l/a = 16, one can 
just about attempt to discern one for l/a = 20 and more easily for higher I. We attempted 
heuristic unstable particle fits to such effective masses, using decay channels weighted with a 
Breit-Wigner for the unstable state as in eqn lfTB"]) . Examples are in Fig. [8] and Fig. [13j The 
interesting feature of such fits is that the true energy is not to be found by looking at E e ff(t) 
at large t, but is typically close to the value of lim^ E e ff(t) i.e. it is very much larger. This 
could significantly reduce the large discrepancy between Nambu-Goto and our observed p = 
excited states. This highlights the dangers in using conventional methods to make energy 
estimates when the state is not stable or very nearly so. 

Despite our uncertainty concerning the existence of many of the excited states, the exis- 
tence of the absolute ground states, whether stable or unstable, appears to be quite unam- 
biguous. And this is often also true of the ground states with p ^ 0. Moreover even when 
the flux is in our rather exotic representations, these particular states have the simple stringy 
excitations of a free string theory, even as / decreases close to its physically minimum value. 
The spectrum of other excited states is more complex and there appears to be plenty of room 
there for the massive modes that might arise from the binding energy in these states. 
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A Representations and Casimirs 



In this Appendix we describe in detail how to calculate flux tubes in the representations of 
interest to us in this paper, and how to determine some of their group theoretic properties. 
Consider a Polyakov loop 

«i =L X 

k= n ( 22 ) 

n x =i 

that winds once around the (spatial) x-torus of length L x in lattice units. (We suppress other 
co-ordinates.) If we take the trace in the representation R, then the operator Tr#{Z p } is a 
candidate operator for projecting onto winding flux tubes carrying flux in the representation 
R. One can deform the Polyakov loop in eqn (l22|) so as to obtain winding operators that can 
be used to form flux tubes with non-trivial quantum numbers. These operators can be used 
as a basis for a variational calculation of the spectrum as described in the text. We use here 
the same basis as we used for the fundamental flux tube in [6], and we refer to that paper for 
a list of the operators. 

The open version of such a flux tube would connect sources in the representations R and 
R. But because the vacuum contains gluons a flux tube in representation R may evolve into 
one in R' by gluon screening if R and R' can be connected by a product of adjoints. And of 
course a flux tube in R may evolve into a product of flux tubes carrying fluxes Ri, R 2 , ... if 
R appears in the product R\ ® R 2 <8> •••• Equally an excited flux tube may decay into a lower 
lying flux tube of the same representation plus colour singlet glueballs. The same remarks 
apply to a closed flux tube winding around the spatial torus. In most cases such a flux tube 
will not be absolutely stable and whether it exists at all, in the sense that an unstable but 
narrow resonance exists, will depend on the dynamics. 

Note that some decays will be suppressed by the usual 'kinematic' large-iV counting ar- 
guments, for example the decay of a sufficiently excited flux tube state into its ground state 
together with a glueball. Other decays are not suppressed in this way, e.g. the decay of 
an adjoint flux tube into a pair of fundamental and anti-fundamental flux tubes. (Although 
it may be that the detailed dynamics suppresses such a decay.) The latter is analogous to 
the expected 'falling apart' as iV — > 00 of a molecular hadronic state (such as, perhaps, the 
inverted nonet scalars at < lGeV), while the former is analogous to the p meson becoming 
stable in that limit. 

Consider the fundamental representation /. Any representation R will appear in a product 
of a number of / and /. Suppose that the number of these is Nf and Nj respectively, then 
we define the A/"-ality of the representation to be 

k = N f - Nj. (23) 

and we use the generic name of a 'fc-string' for such a flux tube. The reason for focusing on the 
./V-ality is that for an SU(iV) gauge theory k is conserved (module N) under gluon screening 
(since gluons are adjoints carrying k — 0). So for N > A the lightest k = 2 flux tube is an 
absolutely stable state, as is the lightest k = 3 flux tube for N > 6. Of course it may be that 
such a k = 2 'flux tube' consists of nothing more than two k = 1 flux tubes. Whether it does 
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or not is a dynamical question which, in fact, has been answered (see e.g. [T6l H]), and we 
know that the lightest k = 2 and k = 3 strings are flux tubes in their own right, which are 
strongly bound, stable states in SU(6) jS]. 

The irreducible representations we consider in this paper include the totally anti-symmetric 
and symmetric of / <8> / (referred to as k = 2 A and k = 2S), and the anti-symmetric, mixed 
and symmetric of / ® / (g> / (referred to as k = 3 A, k = 3M and k = 3S). The construction 
and properties of these are obtained in the Appendices of [16J to which we refer the reader for 
details, in particular for the derivation of the appropriate operators, and for the values of the 
quadratic Casimir, 

C 2 (R) = Ti R T a T a (24) 

where R is the representation and the T a are the generators of the group. 

We also consider in this paper the following representations: the adjoint, which is k = 
and is the non- singlet piece of / <g) /, and the 84 and 120 which are both k — 1 and which 
arise in / <g> / <g> /. These are the representations we shall now describe in more detail. 

A standard and efficient method for dealing with the representations of SU(iV) is Young 
tableaux. For the derivation and rules of use we refer to [25]. We will label a tableau by 
y(Ai,A2, ...) where Aj is the number of boxes in the i'th row of the tableau. We recall the 
rule that Aj > Aj+i, and we only show the \j for rows j containing at least one box. So, for 
example, / corresponds to Y(Xi = 1), / corresponds to Y(Xi = 1, A 2 = 1, Ajv-i = 1), while 
the adjoint is given by F(Ai = 2, A 2 = 1, Ajv-i = 1)- We also recall that the dimension 
of the representation R corresponding to the tableau K(Ai, A 2 , ...) is given by [25] 

d R = * Lh < ' ^— t ; h = X k + N-k , k = l,...,N (25) 

[n — I)' {n — 2J1...1! 

and that its quadratic Casimir is given by [IB] 

-if n r 2 \ 

C 2 (R) = - in b N + J2 Ai(A< + 1 - 2i) - ^ J . (26) 

where = Ai + Ai + ... is the total number of boxes in the tableau and n r is the number of 
rows. 

Applying eqns ()25|l26p to the fundamental tableau, Y(Xi = 1), gives 

df = N ; C7 a (/) = ^^i (27) 

The adjoint representation is obtained from / <g> / = 1 + adj. The singlet corresponds 
to adding the single / box to the bottom of the / column of boxes, while the adjoint is 
obtained by adding it to the first row. Applying eqns( l25|26l) to the resulting adjoint tableau, 
F(Ai = 2, A 2 = 1, Ajv_i = 1), gives 

d ad] =N 2 -l ; C 2 (adj)=N (28) 
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For SU(6) we see that d ac y = 35. 

Consider now the product / <g) / <S> f- As we saw above the first product gives us / <S> f — 
1 + adj. Adding the second / box to the singlet just gives us an / again, while adding it to 
the adjoint Y(Xi = 2,X 2 = l,...,Ajv-i = 1) gives us R a = Y(X± = 3, A2 = 1, ...,Ajv-i — 1) if 
we add it to the first row and R b = Y(Xi = 2, A 2 = 2, , A 3 = 1, Xjy^i = 1) if we add it to 
the second, as well as a second / by adding it to the bottom of the first long column (which 
is then of length N and so can be dropped). Applying eqns( l25|26|) to these tableaux gives, 

«U. = i*(tf + 2)(*-l) ; *<«.) = (J^m±A = C2(/) H_i (29) 

and 

dRt ^ N{N+l){N - 2) ; QW = M±^Zil = ft(/) ^. ,30) 

For SU(6) we can label these representations by their dimensions, i.e. 

R a = 120 ; R b = 84 : for SU(6) (31) 

and evaluating the quadratic Casimirs we obtain the entries in Table EJ 

To calculate the spectrum of a flux tube in representation R we want correlators of op- 
erators Tr R / p where l p is some loop winding (once) around the appropriate spatial torus, the 
simplest example being that in eqn fr2"2~]) . Since our Monte Carlo generates group elements 
Ui that are in the fundamental representation, we need to express Tr R Z p in terms of Trf/ p . 
This can be done by taking products of group elements and imposing the (anti) symmetry 
constraints on the indices that are encoded in the corresponding Young tableaux, as carried 
out explicitly in [16] for the 2A, 2S and 3A, 3M, 35* representations. Often we can employ a 
short cut. For example if we want R and we find that R\ ® R2 = R3 © R then we can use 
the fact that Tr Rl Z x Ttr 2 1 = TrR 3 Z + Tr^l. For example, we know that the adjoint satisfies 
f ® f = adj © 1 (where 1 is the singlet) so 

Tr adj Z = Tr f /Tr F / - 1 = |Tr f /| 2 - 1. (32) 

We also want such an expression for the representations R a and R b . One can easily see that 
2S © / = Ra © / giving 

Tr R J = Tr 2S /Tr f /t - Tr f Z = l - (Tr f 2 Z + Tr f / 2 ) Tr^ - Tr f Z (33) 

and that 2A © J = R b © / giving 

Tr Rb Z = TrsA/Trf/t - Tr f Z = ^ (Tr 2 / - Tr f / 2 ) Tr f /t - Tr f Z (34) 

where we use the fact that Tijl = Trf P and we obtain T^s^a^ from [16] . The above expressions 
are for general iV and for SU(6) give us the appropriate operators for the adjoint, 120 and 84 
respectively. 
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Figure 1: Effective energy of the k=2A, p=0, P=+ (variational) ground state of a flux tube 
of length l/a = 16,20,24,28,32,36,40,44,48,52,64. Lines are our plateaux estimates (±1(7 
error bands). 
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Figure 2: k = 2A ground states with p = 0, 2-n/l, Air /I and with P = +, o, and P = — , •. Solid 
red curves are Nambu-Goto predictions. Dashed red lines are the model in eqn( |T5l) . Dashed 
blue lines denotes lower boundaries of scattering states formed of two fundamental flux tubes 
with total momentum p. Vertical line denotes location of 'deconfinement' transition. 
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Figure 3: Energy of k = 2 A ground state with p = and P = +, minus predictions of various 
'models': Nambu-Goto, •; linear plus Liischer correction, o; and only linear term, □. The 
solid curve includes an 0(1/ 'I 7 ) correction to Nambu-Goto. Vertical line denotes location of 
'deconfmement' transition. 
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Figure 4: Phonon excitation energies, as defined in eqn (|T4|) . of k = 2 A ground states with 
p = 0, 2tc/1, Ait /I and with P = + (•) or P = — (■). Open symbols shown for p = 0,2tt/1 
are without the zero-point energy in eqn fTll"]) . Horizontal lines are Nambu-Goto predictions. 
Vertical line denotes location of 'deconfinement' transition. 
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Figure 5: Fitting the p = 2n/l ground state energies to the model in eqn fll5l) and extracting 
the excitation mass averaged over I > l . For representations r = 2A (•) and r = 3A (o). 
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Figure 6: Energies of lightest k = 2A p = states with P = +, •, and with P = — , o. Solid 
curves are corresponding Nambu-Goto levels. Upper dashed line is (approximately) the energy 
of the lightest P = + decay channel consisting of a flux tube and glueball; lower line is that 
of two fundamental flux loops. Vertical line denotes location of 'deconfinement' transition. 
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Figure 7: Effective energy, E e ff, of the k=2A, p=0, P=+ (variational) first excited state of a 
flux tube of length l/a = 16, 20, 24, 28, 32, 36, 40, 44, 48, 52, 64. Values shown have been shifted 
by multiples of 0.05 for clarity, aE e ff = aE e ff + 0.05 * (l/a — 24)/4, and a shift of 0.35 for 
I = 64a . Lines are ±lcr error bands of our mean plateaux estimates. Points on right axis are 
corresponding decay thresholds (= 2Ef(l)). 
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Figure 8: Effective energy, E e ff, of the k=2A, p=0, P=+ (variational) first excited state 
of a flux tube of length I = 16a. The red curve is what one obtains with an excited state, 
aE* = 0.49, in addition to the ground state, which is indicated by the lower horizontal line. 
The solid black curve is obtained by summing over scattering states of 2 fundamental flux 
tubes with a Breit-Wigner amplitude peaking at the upper horizontal line and with a width 
aT = 0.065. The lower black curve sums uniformly over all scattering states. 
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Figure 9: k = 2S ground states with p = 0, 27r/Z, Att/I and with P = +, o, and P = —, 
• . Curves are Nambu-Goto predictions. Dashed lines are thresholds for scattering states 
formed of two fundamental flux tubes with total momentum p. Vertical line denotes location 
of 'deconfmement' transition. 
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Figure 10: Phonon excitation energies (see eqn( H4p ) of lowest k = 2S states with p = and 
P = +, •, p — 2tt/1 and P — — , o, and with p = 47r/Z for both P = — , o and P = +, •. 
Horizontal lines are Nambu-Goto predictions. Vertical line denotes location of 'deconfmement' 
transition. 
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Figure 11: Effective energy of the k=2S, p=0, P=+ (variational) ground state of a flux tube 
of length l/a = 16, 20, 24, 28, 32, 36, 40, 44, 48, 52, 64. Plateau estimate (±1(7 error band) given 
by red lines. Lowest energy of two fundamental lines is indicated by diamonds on right axis. 
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Figure 12: Effective energy of the k=2S, p=0, P=+ (variational) first excited state of a flux 
tube of length l/a = 16, 20, 24, 28, 32, 36, 40, 44, 48, 52, 64. Lowest energy of two fundamental 
lines is indicated by diamonds on right vertical axis. 
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Figure 13: Effective energy, E e ff, of the k=2S, p=0, P=+ (variational) first excited state of 
a ffux tube of length I = 16a. The red curve is an example of fitting with an excited state, 
aE* = 0.595, in addition to the ground state indicated by the lower horizontal solid line. The 
upper black curve is obtained by summing over scattering states of 2 fundamental flux tubes 
with a Breit-Wigner probability peaking at the upper horizontal line, with width aT = 0.12. 
The lower black curve sums uniformly over all scattering states. 
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Figure 14: k = 3A ground states with p = 0, 2ir/l, An/l and with P — +, o, and P = — , 
• . Red curves are Nambu-Goto predictions. Black dashed line denotes lower boundary of 
scattering state formed of three fundamental flux tubes with total momentum p, and black 
solid line of a k = 2A flux tube with a fundamental flux tube with total momentum p. Vertical 
line denotes location of 'deconfinement' transition. 
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Figure 15: Effective energy of the k=3A, p=0, P=+ (variational) ground state of a flux tube 
of length l/a = 16,20,24,28,32,36,40,44,48,52,64. Lines are our plateaux estimates (±1<t 
error bands). 
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Figure 16: Energies of lightest k = 3 A states with p = and with P — +, •, or P = —, o. 
Red curves are corresponding Nambu-Goto levels. Black dashed line is the lowest energy of 
three fundamental flux loops; solid black line is that of a k = 2A flux tube with a fundamental 
one. Vertical line denotes location of 'deconfinement' transition. 
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Figure 17: Effective energy, aE e ff, of the k=3A, p=0, P=+ (variational) first excited state 
of a flux tube of length l/a = 16, 20, 24, 28, 32, 36, 40, 44, 48, 52, 64. Values shown have been 
shifted for clarity: aE e ff = aE e ff + 0.05 * {l/a — 32)/4. (Shift of I = 64a is 0.25.) Lines are 
our plateau estimates (±lcr error bands). Thresholds for decay shown as diamonds on right 
axis. 
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Figure 18: k = 3M ground states with p = 0, 27r/Z, 47r/Z and with P = +, o, and P = 
— , •. Solid red curves are Nambu-Goto predictions. Dashed line denotes lower boundary 
of scattering state formed of three fundamental flux tubes with total momentum p; black 
line of a k = 2A flux tube with a fundamental flux tube. Vertical line denotes location of 
'deconfmement' transition. 
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Figure 19: Effective energy of the k=3M, p=0, P=+ (variational) ground state of a flux tube 
of length l/a = 16,20,24,28,32,36,40,44,48,52,64. Lines are our plateaux estimates (±1<t 
error bands). Decay thresholds indicated by diamonds on right axis. 
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Figure 20: k = 3S ground states with p = 0, P = +, •, and with p = 2n/l, P = — , o. Solid red 
curves are Nambu-Goto predictions. Dashed line denotes lower boundary of scattering state 
formed of three fundamental flux tubes with same momentum; black line of a k = 2A flux tube 
with a fundamental flux tube. Vertical line denotes location of 'deconfinement' transition. 
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Figure 21: Effective energy of the k=3S, p=0, P=+ (variational) ground state of a flux tube of 
length I /a = 16, 20, 24, 28, 32, 36, 40, 44, 48, 52, 64. Lines are our plateaux estimates (±lcr error 
bands). No plateau attempted for I = 48a. Diamonds on right axis denote decay thresholds. 
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Figure 22: Adjoint ground states with p = 0, 2n/l, 4tt/1 and with P = +, o, and P = — , 
• . Solid red curves are Nambu-Goto predictions. Dashed lines denote lower boundaries 
of scattering state formed of a pair of (anti)fundamental flux tubes with same momentum. 
Vertical line denotes location of 'deconfmement' transition. 
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Figure 23: Phonon excitation energies of adjoint ground states with p = and P — +, •, 
p = 27r/Z and P — — , o, and with p = 47r/Z for both P = —, o and P = +, •. Lines are 
Nambu-Goto predictions. 
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Figure 24: Ground state adjoint loop effective masses for l/a =16, 20, 24, 28, 32, 36, 40, 
44, 48, 52, 64. Solid lines give ±lcr error bands of our plateaux estimates. Dashed lines 
are twice the energy of corresponding (anti)fundamental loops. (Also indicated by diamonds 
on right axis.) Values shown have been shifted by multiples of 0.025 for clarity: aE e ff = 
aE eff + 0.025 * (l/a - 16)/4, except a shift of 0.25 for I = 64a. 
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Figure 25: Ground state adjoint loop effective masses for I /a = 16, 20, 24, 28. Solid lines give 
±lcr error bands for our plateaux estimates. Dashed lines are error bands for twice the energy 
of corresponding (anti)fundamental loops. 
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Figure 26: Ground states in the 84 representation for p = 0,2tt/1,Att/1 and with P = +, 
o, and P = — , •. Solid red curves are Nambu-Goto predictions. Black dashed line denotes 
lower boundary of scattering state formed of three (anti)fundamental flux tubes with same 
momentum, black solid line of one k = 2A and one antifundamental flux tube. Vertical line 
denotes location of 'deconfmement' transition. 
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Figure 27: Ground state effective masses of 84 flux tube with l/a =16, 20, 24, 28, 32, 36, 40, 
44, 48, 52. Solid lines give error bands of our plateaux estimates. Diamonds on axis indicate 
2A + / (solid) and f + f + f (open) decay thresholds. 
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Figure 28: Ground states in the 120 representation for p = 0,2tt/1 and with P = +, o, 
and P — — , •. Solid red curves are Nambu-Goto predictions. Black dashed line denotes 
lower boundary of scattering state formed of three (anti)fundamental flux tubes with same 
momentum, black solid line of one k = 2A and one antifundamental flux tube. Vertical line 
denotes location of 'deconfmement' transition. 
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Figure 29: Overlap, as in eqn(l2ip. of k = 2A and adjoint ground states onto low-lying 
fundamental states. For / = 32a. 
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Figure 30: Overlap, as in eqn(l2ip. of k = 2A and adjoint ground states onto low-lying 
fundamental states. For I = 64a. 
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Figure 31: Overlap, as in eqn( l2ip . of k = 2A first excited p = state onto low-lying funda- 
mental states. For I = 32a and I = 64a. 
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Figure 32: Overlap, as in eqn( l2ip . of k = 3 A first excited p = state onto low- lying funda- 
mental states. For I = 32a and I = 64a. 
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Figure 33: Overlap, as in eqn (!2T|) . of k = 2A second and third excited p = states onto 
low-lying fundamental states. For I = 32a. 
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Figure 34: Overlap, as in eqn (!2T|) . of k = 2 A second and third excited p = states onto 
low-lying fundamental states. For I = 52a. 
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